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ABSTRACT 

An abstract approximation theory and computational methods are developed 
for the determination of opitmal linear-quadratic feedback controls, observers 
and compensators for infinite dimensional discrete-time systems. Particular 
attention is paid to systems whose open-loop dynamics are described by 
semigroups of operators on Hilbert spaces. The approach taken is based upon 
the finite dimensional approximation of the infinite dimensional operator 
Riccati equations which characterize the optimal feedback control and observer 
gains. Theoretical convergence results are presented and discussed. 
Numerical results for an example involving a heat equation with boundary 
control are presented and used to demonstrate the feasibility of our methods. 
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1. Introduction 


In this paper we develop an approximation theory and computational methods for 
the determination of the optimal feedback control law for the discrete-time linear- 
quadratic regulator problem, the optimal state estimator or observer gains and the 
optimal compensator for Infinite dimensional systems. Specifically, we are 
concerned with systems whose dynamics can be described in terms of linear semigroups 
of operators on Hilbert spaces. The essential feature of our approach is the finite 
dimensional approximation of the infinite dimensional operator Riccati equations 
that characterize the optimal feedback control and observer gains. We develop a 
general, abstract approximation framework and an associated convergence theory which 
is applicable to a wide class of problems. 

The theory for the discrete-time control problem has been developed previously 
in [5], while the theory for the observer and compensator in the continuous-time 
case (with particular emphasis on systems describing the vibration of flexible 
structures) is treated in [4]. Along with presenting the theory for the discrete- 
time observer and compensator here for the first time, we briefly review and outline 
our earlier results for the control problem. 

Our treatment below requires that both the discrete-time input and output 
operators be bounded. As will become evident from the example we present in Section 
4 however, an unbounded input operator (i.e. its range is contained in some space 
larger than the underlying state^ space) in the continuous-time problem may lead to a 
bounded input operator when the system is sampled and considered in a discrete-time 
setting. This notion can, to a certain extent, be generalized to permit the 
application of the theory we develop here to a wide class of problems which 
simultaneously involve unbounded input and output operators. This idea will be 
treated in a forthcoming paper. 
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Now we provide a brief outline of the remainder of the paper. In Section 2 we 
present the theory for the optimal infinite dimensional control law, observer and 
compensator. In Section 3 the approximation theory and convergence results are 
discussed. An example (including numerical results) involving a heat equation with 
boundary control is used in Section 4 to demonstrate the feasibility of our methods. 

2. The Infinite Dimensional Optimal Control Law, State Estimator and 
Compensator 

Let {H, <•,•>} be a Hilbert space and consider the time invariant, discrete- 
H 

time linear control system 

(2.1) z k+l = ^ z k + ^ u k* k — 0,1,2,.... 

z o £ H 

(2.2) = Cz^ + Du^, k — 0,1,2,.... 

where T e L (H), Be L(R m ,H), C e L(H,R P ) and D e L(R m ,R p ). The infinite time 
horizon discrete-time linear-quadratic regulator (LQR) problem is given by: 

* * oo 

Find u = {u^J-^q e ) which minimizes the quadratic performance index 

03 

( 2 * 3 ) J (u) = I <Q Z k’ Z k> H + U k Ru k 

k=0 ' *■ 

where Q e L(H) is self-adjoint and nonnegative, R e L(R m ) is a symmetric, positive- 

OO 

definite m x ra matrix, z^ e H is given and z = {z^^o is determined by the 
recurrence (2.1). 
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As in 
frequently 


the finite dimensional case the discrete-time system (2.1), (2.2) is 
the result of sampling a continuous time system of the form 


(2.4) z(t) = Az(t) + Bu(t), t > 0 

z(0) = z Q 


(2.5) y(t) = Cz(t) + Du( t ) , t > 0 

where A is the infinitesimal generator of a Cg-semigroup of bounded linear 
operators, (T(t) : t > 0} , on H and B e L(R m ,H). In this case we have 
T = J(t) and B = Jq T (t)8 dt where x is the length of the sampling interval or 
sampling period. 

A control sequence u £ Z^(0,®;R m ) is said to be admissible for the initial 
conditions Zq if J(u) < 00 . It can be shown (see [5], [10]) that if there is an 
admissible control for each z^ £ H, then there exists a nonnegative self-adjoint 
solution II e L(H) to the operator Riccati algebraic equation 

(2.6) n = t*(h - iib(r + B*nB) -1 B*n)T + q. 


If, in adddition, u admissible for Zq implies lim |z^|^ = 0 then this solution is 

k*°° 

unique. Moreover under the two hypotheses given above, the LQR problem admits a 

k & 

unique solution u for each z n £ H with J(u ) = <IIz ,z n > • The optimal control is 

(J U U H. 

given in feedback form by 


u 


* 

k 



k = 0,1,2,... 


where the optimal feedback control gains F are given by 
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(2.7) 


* -1 * 

F = (R + B nB) B H 


* * 


and z = is the resulting optimal state trajectory. We have 


Z k+l Sz k’ 


z 0 z 0 


k = 0,1,2,. 


with the optimal closed-loop state transition operator S given by 


( 2 . 8 ) 


S = T - BF. 


If Q is also coercive (i.e., Q > a for some ct > 0) then S has spectral radius less 
than one and S is uniformly exponentially stable with 


s k | < (|n|/o)(i - a/ In | ) k . 


k = 0,1,2,... . 


For each j = l,2,...,m an application of the Riesz Representation Theorem 
th th 

yields that the j component of the optimal k L control input is given by 


[uj. = 
K J 


<f j * z k > H 


for some f . e H. The vector f = (f,,f 0 ,...,f ) e x H is 

j m j=i 

optimal functional feedback control gains. 


referred to as the 


In order to implement a feedback control law of the form = -Gz^, k = 
0,1,2,... where G e L (H,R m ) it Is necessary that the full infinite dimensional state 
z, be available for each k. In practice, however, only a finite dimensional 
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observation e R P of the state, as given in (2.2), is provided. Consequently, a 
state estimator or observer is required. 

A 

For any operator G e L(R P ,H) the discrete-time linear system 


(2.9) 


'k+1 


Tz k + Bu k + G{y k " Cz k - V* 


s H 


k = 0,1,2,... 


is called an observer or estimator for Che system (2.1), (2.2). The feedback 
control law 


A 

(2.10) u k = -Gz k , k = 0,1,2,.. 

along with the observer (2.9) is referred to as a compensator for the system (2.1), 

( 2 . 2 ). 

A 

If we define e fc = z fe - z k , k = 0,1,2,... then direct calculation yields 

A 

e k+l ~ ^ GC) e k> k = 0,1,2,... 

^ ^ k ^ A 

or = S(G) , k ■ 0,1,2,... where S(G) e L(H) is given by 

A A A 

(2.11) S(G) - T - GC. 


If the control law or compensator (2.10) is to be used, then it is desirable to 

A 

have 0 as k ■* 00 . The observer corresponding to G is said to be strongly 

stable if |S(G) Zg|^‘^ 0 as k> 00 for each z^ e H. It Is said to be uniformly 

A A A 

exponentially stable If there exist positive constants M and r with r < 1 such that 


If z = { z k > k _ 0 and z = { z fc } k=Q are generated by (2.1) and (2.9) respectively 

oo 00 A X 00 

with u = {u k ) k= o given by (2.10) then z = tz k ^ k= Q = ^ z k ,Z k^ ^ k=0 satisfies the 

recurrence 

(2.12) z k+l = S(G ’ G) Z k’ k= 0,1,2,... 

A 

where S(G,G) £ L(H x H) is given by 

T 

S (G,G) = A 

GC 

The system (2.12) or equivalently 

z k = S(G,G) k Zq, k- 0,1,2,.. 

is the closed-loop system corresponding to the control system (2.1), (2.2), the 
observer (2.9) and the compensator (2.10). 

We recall (2.8) and by analogy to (2.11), for G e L(H,R m ) we adopt the notation 

(2.13) S(G) = T - BG. 



Using the facts that ' 
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z k+l ~ Tz k BGz k 

= (T - BG)z k + BGe k 

= S(G)z k + BGe k 

and 

e k+1 = [I, -I] z k+1 = [I, -I] S(G,G) z k 
" S(G) [I , -I] z k 
= S(G)e k 


and consequently that 


S(G,G) = U 


S(G) 

0 


BG 

S(G) 


U 


-1 


where U = U * e L(H x H) is given by 


U = 


0 

-I 


standard arguments can be used to establish the following result. 


Theorem Suppose that there exist positive constants M, M, r and r for which 


| S(G; | < Mr and- |S(G) | < Mr 


k = 0,1,2,. 


Then for each n with a. > max(r,r) there exists a positive constant M such that 

•|S(G,G) k |< M* k , 


k = 0,1,2,*.. 
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In particular if S(G) and S(G) are uniformly exponentially stable (i.e. 

A A A A 

r, r < 1) then so too is S(G,G). Also, the spectrum of S(G,G), o(S(G,G)) is given 

by 


a(S (G,G)) = a <S (G) ) U a(S(G)). 


By analogy to the finite dimensional case (see [7]) we define the optimal 
discrete-time observer for the system (2.1), (2.2) to be the system (2.9) with the 

A A 

observer gains G replaced by F given by 


F = me (R + 


A * -1 

cnc ) 


where II e L (H) is the minimal nonnegative self-adjoint solution (if one exists) to 
the Riccati algebraic equation 


(2.14) 


A A A ^ A A ^ __ I A ^ A 

n = T(n - nc (r + cue ) oi)t + q , 


Q e L (H) is nonnegative self-adjoint and R e L (R^) is a symmetric positive 
definite p*p matrix. When G in (2.10) is taken to be the optimal feedback control 

A A 03 A * A *00 

gains F given in (2.7) and z = { z^} is taken to be z = > t * ie trajectory 

determined by the optimal observer, the resulting feedback control law 


(2.15) 


u fc = -Fz fc , k = 0,1,2,.. 


is known as the optimal infinite dimensional compensator. The optimal closed-loop 
system is given by 
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Z k+1 = SZ k " S(F ’ F) Z l ’ k = °» 1 » 2 >“ ‘ 

\ 

We note that the adjoint of the optimal observer gains F are the optimal 

I 

f control gains for the linear regulator problem obtained by replacing the operators T 

* * 

and B in (2.1) with the operators T and C and the operators Q and R in (2.3) with 

A A 

the operators Q and R. Consequently the necessary and sufficient conditions for the 

[ 

j existence of a nonnegative self-adjoint solution (and therefore a minimal 

nonnegative self-adjoint solution) H to the Riccati algebraic equation (2.14) become 
clear and can be found in [5]. In addition, it is now also easy to specify the 
conditions under which 1) (2.14) has a unique nonnegative self-adjoint solution and 
2) the operator 


S = S(F) = T - FC 


will be uniformly exponentially stable. 

The optimal observer gains F Is an element in L(R^,H). They therefore have a 
representation of the form 


where f^ e H, i 


Fy = L 

i=l 


1,2,... ,p. 


A T D 

f i y i ’ y = ^ y l» y2’*“ ,3 V E RP 

a ^ A A rjr 

The vector f = (f ^ ,f^ . ,f^) e 


P 

x H is referred to as 
1 


the optimal functional observer gains. 

Ordinarily of course, the otpimal observer has a stochastic interpretation. 
The optimal observer as we have defined it above is the natural extension to 
infinite dimensions of the well known finite dimensional discrete-time Kalman-Bucy 
filter for the case in which the state and output equations are corrupted by 
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uncorrelated, stationary Gaussian white noise processes with zero mean and 

A A 

respective covariance operators (matrices) Q and R (see [7]). Under appropriate 
additional hypotheses (i.e. that the state weighting operator Q in (2.3) is trace 
class, see [1]) it is also possible to provide the optimal infinite dimensional 
compensator given by (2.15) with the standard finite dimensional stochastic 
interpretation (i.e. as the usual optimal LQG compensator, see [3], [4]). The ideas 
which have been presented above and the approximation theory which will be described 
in the next section require only that the conditions which have been set forth thus 
far hold. We shall therefore continue to take a strictly deterministic approach and 

A /K 

assume that the operators Q and R (as well as Q and R) are determined by engineering 
criteria (for example, stability margins, robustness of the closed-loop systems, 
etc.) rather than via an assumed noise model incorporated into the underlying 
dynamics. 

3. Approximation and Convergence 

In this section we develop an approximation framework which yields finite 
dimensional approximations to the optimal infinite dimensional control and observer 
gains and the optimal infinite dimensional compensator. Central to our approach are 
finite dimensional approximations to the infinite dimensional operator Riccati 
equations (2.6) and (2.14). The approximating equations can be solved using 
conventional techniques (for example, eigenvalue-eigenvector or Schur decomposition 
of the associated Hamiltonian matrices). Under mild and rather general assumptions, 
the convergence of the approximation can be argued. 

For each N = 1,2,... let be a finite dimensional subspace of H. Let 
: H * denote the orthogonal projection of H onto with respect to the H- 
inner product, <*,*> H • For s L (H^) , e L(R m ,H N > and Q N e LO^) with Q N 
nonnegative self-adjoint, we consider the Riccati algebraic equation 
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(3.0 n „ - <(!I S - V n <R ♦ +Q(j 


We assume that for each N, the equation (3.1) has a unique nonnegative self-adjoint 
solution n N e LCH^) and define the approximating optimal control gains by 


F » ■ (R + W. 1 ' 1 ™ 


For the estimator, we take e L (H^) nonnegative self-adjoint and 

C T e L(H,R P ) and consider the equation 
N N 


( 3 * 2 ) n N T N (n N - n N c N (R + <yi N c N ) + Qn 


Assuming a unique nonnegative self-adjoint solution e 1(11^), we define the N" 
approximating optimal observer gains by 


F N T N II N C N ( ' R + C N II N C N^ 


The N l approximating optimal compensator is given by 


(3.3) 


U N,k -F N Z N,k’ 


k = 0,1,2,. 


where z^ = {z^ ^ determined from the N approximating optimal observer 


( ' 3 ' 4) Z N,k+l T N Z N,k + B N u N,k + F N {y N,k °N Z N,k ^N.k* ’ k °» 1 > 2 ’' 


'N,0 P N Z 0 £ ^ * 
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* OP 


The output sequence y N = { y N k > is given by 


J l,k ~ Cz N,k + ““B.k 


k - 0,1,2,. 


where 


Z N,k+l Tz N,k " BF N Z N,k * 


k = 0,1,2,... 


"n.o z o * 


The approximating optimal closed-loop sytem evolves according to the recurrence 


* * 

z = S 


N, k+1 ” S N ^N, k’ k " 


" here 4,k ■ (z »,k- Z B,k )T and S » c L<H * V 18 glve " by 


5 b' 


V 


-BF. 


t n - b n f n - f n c 


N 


The equations and formulas given above are operator equations and as such can 

not be used in computations directly. It is their matrix equivalents with respect 

to a given basis for H„ that are required. 

N ^ 

We assume that the collection {<J> is a (not necessarily orthogonal) basis 

and define $ e * H 

1 

notational convention that for a linear operator L with domain and range in Hjj, R m 


for H n 


. * N /jN . N 
N by $ * C+j, + 2 ,... 


<Ptr ) • 


We adopt the 
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d N N 

or R H , its matrix representation with respect to the basis for H N and the 

standard bases for R ra and RP will be denoted by [L ] . 

If we define the Gram matrix vF = <$^,($^)^>^ then [T^ ] = (M N ) ^[T^] T M N , 

[B* ] = and [C* ] = (M N ) _1 [C n ] T . Also r N = M N [II N ] and r N = (n N ](M N ) _1 

are respectively, the unique nonnegative symmetric solutions to the matrix Riccati 

algebraic equations 

N TNN TN-lTN NT 

(3.5) r w = [t n ] (r - r N [B N ](R + [b n 1 J r w [ b n ] ) [b n ] 1 r N )[T N ] + q n , 


and 


(3.6) r N = [t n ](? n - r N [c N ] T (R + Cc N ]r N [c N ] T ) _1 [c N ]r N )[T] T + q n 


where = M^[Q^] and = [Q^](M^) ^ . The matrix representation for the 
approximating optimal control gains is given by 


[f n ] = (R + [B N ] T r N [B N ]) l [B N ] T r N [T N ] 


and for the approximating optimal observer gains by 


[F N ] = [T N U N [C N ] T (R+ [c N ]f N [c N ] T ) _1 . 


A tAr m ip a ^ A 

If we write ^ ^ with ^ e R ,k = 0,1,2,... then from (3.3) we 


U N,k _[ V ? N,k’ 


k = 0,1,2,. 


obtain 
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and from (3.4) 


S,k+1 ^ T N^N,k + ^ B N^ U N,k + ^ F N^ y N,k“ *- C N ^ C N,k ’ k O.i.Z,... 


A * /W N .-1 N 
C N,0 ^ <$ ’^H * 


m 


If we let f N = (f^, f»,...,f N ) T e * H. T denote the N tk approximating optimal 
I Z m ^ N 

functional feedback control gains, then from 


F N Z N <f ’^H ^ F N^N 


N T 

where z._ = ($ ) r e H T , we find 

N N N 


-N rT , 1A/ N r LN 

f - [f n ](m ) $ 


Similarly, the N L approximating optimal functional observer gains 

;n ,;n :n :n,t , 

f = (f 1 ,f 2 ,... ,f p ) e x Hjj are given by 


f N = [f n ] t $ n 


It is immediately clear that the limiting behavior of the approximation (i.e. 
as N + ») is determined by the limiting behavior of the solutions to the finite 
dimensional Riccati equations (3.1) and (3.2). A convergence theory for 
approximations to discrete-time Riccati equations was developed in detail in [5]. 
We briefly summarize those results here. 

We assume that the spaces H N are H-approxi mating in the sense that the 
projections P N converge strongly to the identity on H as N + 00 . Also, we require 

£ £ A A 

that T n P n z -*• Tz, T^z + T z, Q^z ♦ Qz and + Qz for each z e H 
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and that B ■> B and C P > C in norm as N + ® • Define S XT , S._ e L (H XT ) by 
N NN N NN 

S N = T N- B N F N and S N = T N " F N C N * 


If there exists a positive constant M (M) independent of N for which 


11^ < M (II N < M) then there exists a nonnegative self-adjoint solution II (II) of (2.6) 

A A 

((2.4)) and (II converges weakly to II (II) as N + » . If, in addition, there 

A 

exists a positive constant r (r) less than one and independent of N for which 


(3.7) 


< M r 


< M r ), 


k » 1,2,, 


then n N P N will converge strongly to II (II). 

A 

If the operators (Q^) are uniformly (with repsect to N) coercive and 

A a 

the II N (11^) are uniformly bounded then there exists a positive r (r) less than one 

A 

for which (3.7) holds. If it is also true that Q (Q) is trace class and 

A A A 

^N P N conver g es * n trace norm to Q (Q), then II (II) is also trace class and 

a A 

H N P N (n N P N ) converges in trace norm to H (II). 

The consequences of these results in the context of the control and observer 

problems are at once clear. If H X P _ H weakly as N 00 , then 

N N 

N 

■* F and S^P^ * S strongly and f^ -► f ^ , I = weakly in H as 

N ■* °°. If II^P^ * H strongly, then F^ + F in norm, S^P^ + S strongly and 
N 

f^ + f ^ , i = l,2,...,m strongly in H as N > °°. For the observer problem, 

A A AAA A A^A 

if n N P N + n weakly, then F N F and S N P N S weakly and f -*■ f , i = l,2,...,p 

A A A A A A 

weakly in H as N -*• If H P -*■ n strongly, then F„ F in norm, S. T P. T S 

N N N N N 

strongly and f^ f^, i = l,2,...,p strongly in H as N 00 . 

Let P N denote the projection of H x H onto H x defined by 
P N (z 1 ,z 2 ) = ( z i» p n z 2^* If n N p N n weakly or strongly then 

A A 

+ S weakly or strongly depending only upon whether H N P N n weakly or strongly 
as N + “ . Under appropriate additional hypotheses on the spectral properties of 
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the open loop system and the nature of the approximation spaces and the 
projections P^, it is possible to obtain a result regarding the norm convergence 
of to S (see [4]). An important consequence of this norm convergence of the 

closed-loop systems is that the uniform exponential stability of S would imply the 
uniform exponential stability of for all N sufficiently large. 

* * 

Remark It is sometimes the case that while T^P^ *► T strongly, T^P^ * T 
only weakly (see, for example [2]). In this instance it remains possible to 
demonstrate the weak convergence of II to n (see [5]). However, we do not see how 
to retain the weak convergence of A ^P^ to It. 

Remark If the discrete-time system (2.1), (2.2) was obtained via the sampling of a 
continuous time system of the form (2.4), (2.5) then the approximation to T, T^, is 

frequently obtained by approximating the operator A by an operator A^ and then 

* * 

setting = exp (A^t). The convergence of to T and T^P^ to T can then be 

argued using the well known Trotter-Kato semigroup approximation result (see [6]). 

4. An Example: A Heat Equation with Boundary Input 

We consider the parabolic system with boundary control given by 

Tt = §7 a ( x ) t > 0, X e (0,1) 

w(t,0) = 0 w(t,l) = v(t), t > 0 

w(0,x) = <p (x) , x e [0, 1 ] 

where a e H 1 (0,1), a(x) > 0, x e [0,1], <J> e L 2 (0,1) and v e L 2 (0,~). We take the 
average temperature over an interval of small, but positive, length. 
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[e i ,e 2 ] - That is 

1 ^ o 

y(t) = _ / w(t,x)dx, 0 < e < e_ < 1. 

e 2 £ 1 £ 1 1 

We choose the state space H to be L 2 (0,1) endowed with the usual inner product 
and denote the length of the sampling interval by x. We consider piecewise constant 
controls of the form 

v(t) = u^, t e [kr,(k+l)x), k = 0,1,2,... 

and take the discrete-time state z^ e L^COjl) to be 

z = lim _ w(t,*), k = 1,2,... 

k t+kr 

z 0 = + • 

The resulting discrete-time control system is given by (see [5]) 

z k+1 = Tz^ + Bu^, k — 0,1,2,... 

z 0 = * 

y k = Cz k» k = O. 1 . 2 *”* • 

The open-loop state transition operator T is given by T = T (x ) where 

(T(t) : t > 0} is the analytic semigroup of bounded linear operators on H with 

infinitesimal generator A defined by = (a»|> ' ) ' for 

<|> e H 2 (0,1) fl Hq ( 0 , 1 ). The input operator B e L(R,H) = H is given by 



B - ((I - T(t))* 0 + f X Q T (o )a 'da 


where 4 *q e H is given by 'I'q(x) = x, x £ [0,1]. The output operator 
C e L(H,R) = H' takes the form 

1 ^9 

Cty = -- -- - - /_ *Kx)dx, <J» e L 9 (0,1). 

£ 2 £ 1 £ 1 

The performance index for the control problem is assumed to be of the form 

00 

2 2 

J(u) = l q|z | + ru 

k=0 

A 

with q > 0 and r > 0. For the optimal observer problem we assume that q > 0 

A 

and r > 0 are given. 

The operator A is densely defined, self-adjoint and has compact resolvent, 
satisfies the dissipative inequality 

(4.i) <am> h < -w |* ||[ , * e h 2 (o,d n hJ(o,d 

for some ca > 0 and consequently the semigroup {T(t) : t > 0} is uniformly 
exponentially stable with |T(t)| < e"^, t > 0. 

It follows therefore, that both the optimal control and observer problems 
(along with the associated operator Riccati equations) have unique solutions, 
optimal control law is given by 

u v = “ Fz v = ■ <f » z t > H = ~ In f(x)z*(x)dx, k = 0,1,2,... 
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where f e L^COjl) is the optimal functional feedback control gain. The optimal 
observer gains have the form 


Fy = fy, y e R 

A 

where f e Is the optimal functional observer gain. 

We note that if a(x) » a, a constant, and e ^ and are chosen appropriately, 

then all of the open-loop modes will be controllable and observable (i.e. B and C 

are not orthogonal to any of the eigenfunctions of T = T ). 

We use a standard RItz-Galer kin approach to define a linear spline based 

N N— 1 

approximation scheme. For each N - 2,3,... let = 1 denote the usual ”hat M 

functions on [0,1] which vanish on the boundary. They are given by 


Nx - j + 1 
4> j (x) = \ j + 1 ~ Nx 


z± 

N ’ N 

■JL i +1 


x e [*VS 


N 


xe Ip 

elsewhere 


for j = 1,2,... ,N-1 . 

Letting V be the space Hq( 0,1) endowed with the inner product 

<^l,^2>v = usua -^ com P act embeddings V c H c V* hold. Set 

= span 4 N } N ~! and denote by P N the orthogonal projection of H onto H N with 

respect to the H— inner product. Note that Hjj c V and denote by the orthogonal 

projection of V onto with respect to the V inner product. 

From (4.1) we find 0 e p(A). Consequently A * exists and is compact. 

—1 V —1 1 

Define A„ e L(H„) to be the inverse of the operator A = P M A . It is not 

N N N N tL, 

difficult to show (see [5]) that the operator A N is well defined, self-adjoint and 


satisfies 
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< Vn-+n > h ‘ ■“I*hIh> *N eH H- 

Setting T^(t) - exp (A^t), the semigroups of bounded linear operators on 
H n , {T^(t) : t > 0} are uniformly exponentially stable with ( T^( t ) | < e" Qit > t > 0, 
Elementary properties of linear spline functions (see [9]) imply that 
P -► I strongly on H and > I strongly on V as N + °°. Since A ^ is compact 
and 


|P l A ^ - A ~^I h < l P N A "^ ' A ^Iv = l (P N “ I)A ~^! V ’ * £ H 


V i-1 

we conclude that P* T A 
N 


A 1 in norm as N + 00 . It follows that A„^P„ + A ^ 

N N 


strongly as N > 00 and therefore (using the Trotter-Kato semigroup 

* • k 

approximation theorem, see [6]) that T (t)P XT j> -*• T(t)ty and T x ,(t)P,J> T (t)>|> 

N N jN N 

as N > « for each ip e H uniformly on bounded t-intervals. 


Setting C xt CP m an d 


N y H N 


' N* N 


b n . a-v)^ 0 *r ) T,wp ia '», 


the uniform exponential stability of T N implies that the N th approximating 
Riccati equations (3.1) and (3.2) have unique, nonnegative, self-adjoint 

A A A 

solutions II and II for each N and that IL 7 P^ *► H and H P._ II strongly as 
N N NN NN 

N * 00 • For ^ e H and y e R we have that the approximating optimal 

feedback control and observer gains F^ and satisfy 

V N * = <f N ,P N ^> H = - /J f N (x)^(x)dx 

and 
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for some f N , f N e with -*■ f and -*• f strongly in I^O.l) as N 00 . 

N N 

Recalling the definition of $ and M , the (N-l) x (N-l) matrix 

representation for the operator A^ is given by [A N ] = (M N )~*(L^) where 

L N = - <$ N ,($ N ) T > V . Then [T N ] = exp ( [A N ] t) and defining 

, a' = and 1^ to be the (N-l) x (N-l) identity matrix 

U U rt N rl 

we have 

[B n ] - (I N - [T n ])(M N ) _ Sq + /q ex P ([A N ]a)(M N ) - 1 a;da 

- (I N - [T N ])(M N r 1 ^ + [A n ] _ 1 ([T n ] - I N )(M N ) _ 1 a; , 

[Q„] = ql , [Q n ] = ql and [C N ] = C(* ) with f = [F N ] (M ) V and 

JN rv i t a n 

f = [f n ] * . 

A <*, 

Setting a(x) = 1, x e [0,1], q = 1, q = 1, r * 1, r = 1, x = .01, 

= y - .04/2 and z ^ = y + .03/2 we used our scheme to obtain the 

N A N 

approximating functional gains f and f for various values of N plotted in 
Figures 4.1 and 4.2 respectively below. The matrix Riccati equations (3.5) 
and (3.6) were solved using a generalized eigenvector approach (see [8]). All 
computations were carried out on v an IBM PC personal computer. 
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